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We consider a statistical mixture of two identical harmonic oscillators which is characterized by 
four parameters, namely, the concentrations (x and y) of diagonal and nondiagonal bipartite states, 
and their associated thermal-like noises (T/a and T, respectively). The fully random mixture of two 
spins 1/2 as well as the Einstein- Podolsky- Rosen (EPR) state are recovered as particular instances. 
By using the conditional nonextensive entropy as introduced by Abe and Rajagopal, we calculate the 
separable-entangled frontier. Although this procedure is known to provide a necessary but in general 
not sufficient condition for separability, it does recover, in the particular case x = T = (Va), the 
1/3 exact result known as Peres' criterion. This is an indication of reliability of the calculation of 
the frontier in the entire parameter space. The x — frontier remarkably resembles to the critical 
line associated with standard diluted ferromagnetism where the entangled region corresponds to the 
ordered one and the separable region to the paramagnetic one. The entangled region generically 
shrinks for increasing T or increasing a. 
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Quantum entanglement is a quite amazing physical phenomenon, and has attracted intensive interest in recent 
years due to its possible applications in quantum computation, teleportation and cryptography, as well as to its 
connections to quantum chaos flHl9|. A nonextensive statistical mechanics was proposed in 1988 by one of us, 
and is currently applied [pT|— p6| to a variety of thermodynamically anomalous systems which, in one way or another, 
exhibit (multi)fractals aspects. Among these anomalous systems, a prominent place is occupied by systems including 
long-range interactions and Levy distributions. So being, it is after all not surprising that this thermostatistical 
formalism has interesting implications p7|-p9t in the area of quantum entanglement and its intrinsic nonlocality, thus 
showing the confluence of two concepts coming from distinct physical areas. 

Quantum systems can be more or less entangled, which makes relevant the discussion of whether a given system 
is or not separable. Separability, which we shall define in detail later on, is a crucial feature in the discussion on 
whether a quantum physical system is susceptible of a local realistic description with hidden variables. These issues 
were first discussed in 1935 by Einstein, Podolsky and Rosen (EPR) |3(| and by Schroedinger fl3l| , and since then by 
many others JT| — p~9|] . As we have mentioned above we shall see that entropic nonextensivity provides a path through 
which it is possible to discuss quantum entanglement [p7| , p9||3^1 . 

In this work we analyze a composite system of two harmonic oscillators with identical energy spectra. Following 
along the lines of Abe and Rajagopal |27|, we use the nonextensive entropy S q = 1 ~ T ' 1 P (q 6 72.; S\= — Tr pin p), p 



being the density matrix, in order to study the frontier between separable {quantum nonentangled) and nonseparable 
(quantum entangled) regions. More precisely, we determine a frontier in some parameter-space which is either the exact 
one or an overestimation of the separable region. We remind that the entropic arguments that are used within this 
approach provide, like the Peres' partial transpose method [ fL2|| , necessary but not sufficient conditions for separability. 
The mixed state we shall consider for this bipartite system involves Boltzmann-like probabilities in such a way that 
the bipartite spin 1/2 system is recovered as a particular case. More precisely, the 1/3 Peres' criterion for separability 
will emerge as the T = limit of this specific mixed state. An increase of the temperature T incorporates higher 
energy levels in the mixed state under consideration, which has the effect of enlarging the separable region. 

In order to make this work self-contained we will present, in the following, the connection between the nonextensive 
statistical mechanics and quantum entanglement. Let us begin with the definition of a generalized entropy for a 
quantum system 

S q = l ~ T _ r f (<1 G »; Trp = 1; Sx = ~Tr p\np) , (1) 

where p is the density operator. Let us assume now that p is the density operator associated with a composed system 
A + B. Then the marginal density operators are given by pa = Ttb p and ps = TrA p (TrA pa = Trs Pb = 1). The 
systems A and B are said to be uncorrelated (or independent) if and only if 
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p = p A ® PB ■ (2) 

Otherwise they are said to be correlated. Two correlated systems are said to be separable (or unentangled) if and only 
if it is possible to write p as follows: 

w w 
P = Y,P*Pa®Pb (Pi>0Vi;J2Pi = 1 )- (3) 

i=l i=l 

The limiting case of independency is recovered for certainty, i.e., if all {pi} vanish excepting one which equals unity. 
Nonseparability (or entanglement) is at the basis of the amazing phenomena mentioned before and, as already pointed, 
at the center of the admissibility of a local realistic description of the system in terms of hidden variables. As it is 
known correlation is a concept which is meaningful both classically and quantically. Entanglement, in the present 
sense, is meaningful only within quantum mechanics. The characterization of quantum entanglement is not necessarily 
simple to implement, since it might be relatively easy in a specific case to exhibit the form of Eq. (3), but it can be 
nontrivial to prove that it cannot be presented in that form. Consequently, along the years appreciable effort has been 
dedicated to the establishment of general operational criteria, preferentially in the form of necessary and sufficient 
conditions whenever possible. Peres Jl2| pointed out a few years ago a necessary condition for separability, namely 
the nonncgativity of the partial transpose of the density matrix. In some simple situations (like the simple mixed 
state of two spin 1/2) Peres' criterion is now known to also be a sufficient condition. But, as soon as the case is 
slightly more complex (e.g., 3 x 3 or 2 x 4 matrices) it is known now to be insufficient. Within this scenario Abe 
and Rajagopal |27j] recently proposed a different condition, claimed to be a necessary one, based on the nonextensive 
cntropic form S q given in Eq. (1). Let us summarize the idea. The quantum version of conditional probabilities is 
not easy to formulate in spite of being so simple within a classical framework. The difficulties come from the fact 
that generically p does not commute with either pa or ps- Consistently, the concept of quantum conditional entropy 
is a sloppy one. Abe and Rajagopal suggested a manner of shortcutting this difficulty, namely through the adoption, 
of the following definitions of the conditional entropies S q (A\B) and S q (B\A), respectively given by 

S « iAlB) = 1 + (1- 9 )S,(B) ' (4) 

and 

g (b\ a\ — Sg(A + B) - S q {A) 

Sq(BlA) = l + {l-q)S q {A) ' () 

where S g (A + B) = S q (p), S q (A) = S q (pA) and S q {B) = S q (pB). Obviously, for the case of independence, i.e., when 
P = PA® Pb, these expressions lead to S q (A\B) — S q (A) and S q (B\A) = S q (B), known to be true also in quantum 
mechanics. 

Both classical and quantum entropies S q (A + B), S q (A) and S q (B) are always nonnegative. This is not the case of 
the conditional entropies S q (A\B) and S q (B\A), which are always nonnegative classically, but which can be negative 
quantically (see also p3|]). It is therefore natural to expect that separability implies nonnegativity of the conditional 
entropies for all q. This is the criterion proposed in [27|| . 

Peres' criterion seems to be less restrictive than the entropic one in the sense that it might provide a larger quantum 
entangled region. Although we have no proof that it cannot be the other way around, we have not encountered such 
an example. In several cases, including some of increasingly many spins, both criteria have produced the same result, 
presumably the asymptotically exact answer. Details can be found in |Jl^-|l|,|l|,|| . As a general trend, it seems that 
the transposed matrix criterion overestimates the separable region (i.e., underestimates the quantum entangled region) 
less or equally than the conditional entropy criterion does. It could well be that whenever both criteria produce the 
same result, this result is the exact answer |p5| . 

Equations (^) and (j^) have in fact been proved for a classical system fl32f , not for a quantum one. It seems however 
plausible [^7f that they preserve the same form in both classical and quantum cases. Such conditional entropies enable 
what is a necessary condition for separability. Imposing the conditions S q (A\B), S q {B\A) > for all values of q, the 
particular space of parameters becomes divided into two regions: one where any conditional entropy S q is positive 
and the other one - a domain of entangled states - where at least one conditional entropy is negative for some q. 
Therefore, a line (or a frontier, generally speaking) like a critical one, emerges. 

Within this standpoint, it becomes clear that the Boltzmann-Gibbs-Shannon entropy (q = 1) is a concept too poor 
for properly discussing quantum entanglement, a conclusion recently reached also by Brukner and Zeilinger p8fl from 
a different path. An appreciable amount of arguments are now available in the literature which connect quantum 
entanglement and thermodynamics P,p^pl]Jl^-[l7|,p8| . 
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In what follows we consider, for the two oscillators A and B, the basis j^^lm) b = \n, m), where n, m = 0, 1, 2, .... 
Let us define the symmetric and antisymmetric states \n,m ± ) = -^(\n,m) ± \m,n)) for n > m. These states will 
play a role similar to that of the EPR state. The mixed state we are interested in is the following one: 

l-x-y 



PA+B = 



|0,0)<0,0| + |1,1)<1,1| + |1,0)(1,0| + |0,1)(0,1| 

oc n— 1 



+x aJ2e- 2na/T \n,n)(n,n\ + y bJ2 ^ {n+m)/T \n, mT) (n, m~ | , 



(6) 



n=2 



n— 1 m— 



where T > 0, a > 0, 



a ee l/]T e - 2 " Q / T = (1 - e- 2Q/T ) e 4Q 



IT 



71=2 



and 



oo n — 1 



-(n+m)/T 



= 2(1 - e- 1/T )smh(l/T) 



(7) 



(8) 



n— 1 m— 



We easily verify that TrpA+B = 1- Since the eigenvalues of p^+s must be numbers within [0,1], we have that 
< x,y < 1 and x + y < 1. Parameter a measures the importance of the noise associated with the nondiagonal terms 
with respect to the noise associated with the diagonal ones. Parameter T is a temperature-like measure of noise. 
Notice that, for x = y = 0, we have 

1 



PA+B 



|0,0)(0,0| + |1,1)(1,1| + |1,0)(1,0| + |0,1}(0,1| 



1 . 

= 4 U < 



(9) 



where I4 is the 4-dimensional identity matrix. In other words, this state is isomorphic to a fully random state of two 

spins 1/2. Notice also that, for T = and (x,y) = (0,1), we have p A +B = |1, 0~ > (1, 0~ | = 4j |1, 0)(1, 0| - |0, 1)(0, 1 

which is isomorphic to the celebrated EPR state. The separable-entangled frontier for the T = x = particular case is 
known [||j27|j2|,||, namely y = 1/3 (Va). Before we address the discussion of the full frontier in the (x, y, T, a)-space, 
it is convenient to rewrite Eq. m) in the following equivalent form: 



PA+B 



l-x-y 



4 

OO 



[|0,0)(0,0| + |1,1}(1,1| + |1,0 + }(1,0+| 
1 — x — y 



x a^V 2nQ / T |n,n)(n,n| + 
00 n— 1 

y bJ2 E z~ (n+m)/T \n,m-){n,m- 

n—2 m— 



+ y6e- 1 / r )|l,0-)<l,0-| 



(10) 



Due to the orthonormality of all the states appearing in this expression, the calculation of (pA+B) q is easily carried 
out. It follows 



l — x — y\i x q (1 



-2a/T\q 



I _ e -2qa/T 



l-x-y 

4 

(1-- 



•2y(l 



-1/T\ P -1/T 



sinh(l/T) 



+ 



(2y (1 - e- 1 '?) e-^ T sinh(l/T)) (1 + e 



l/T _ e -2g/T) 



2sinh(g/T) (1 - e~il T ) 
Let us now address the calculation of (pa) 9 - We have 

l ~ x - y '. ye - l ' T smh{llT))\Q){0\ 



(11) 



Pb = Tr A p A +B = 



l-x-y 



y sinh(l/T) (e 



-l/T_ e -2/T + e -3/T- 



|1>(1| 



(x (1 - e- 2Q / T ) e - 2 (»- 2 W T + y sinh(l/T) ( e -"/ T - e - 2 "/ T + e -( 2n + 1 )/ r )) |n)(n| , (12) 
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where |n) denotes the n state of oscillator B. It now follows 



Tr( PB ) q = ( 1 \ V - +ye-V T S mh(l/T)V 
1-x-y 



+ { V + y sinh(l/T) {e- x ' T - e- 2 ' T + e~ 3 / T )) 

OO 

+ e- 2a/T ) e- 2 ^ a / T + y sinh(l/T) [eT n ' T - e- 2n ' T + e -(2»+D/T)) * (13) 

n=2 

Now that we have Tr(pA+B) q and Tr(ps) q as explicit functions of (x, y, T, a), we can immediately obtain S q (A\B)(— 
S q (B\A)) through the use of Eqs. (Q) and (§) yielding 

<S g (A|B) as a function of (x, y) for typical values of {q,T) and a = 1 is represented in Figs. 1 and 2. By comparing 
these figures, notice the nonuniform convergence at T = 1/q = 0. 

Fig. 3 exhibits, in (y, g) space, the line verifying S q (A\B) = 0, for particular values of (x, T, a). The physical region 
is the one satisfying < x + y < 1. Given x, for y above this critical line and within the physical region y < 1 — x, one 
has ^(^IB) < 0, which signals the existence of an entangled state. Notice that for sufficiently small x, the strongest 
restriction to the separable region is imposed at q — > oo. As x increases, the minimum of y occurs at a finite q, q m in 
(the possible nonmonotonicity of S q {A\B) with regard to q has already been noticed by p^j). Since the strongest 
restriction to separability corresponds to the nonnegativity of both S qmin {A\B) and S qmirl (B\A), then this condition 
is the one which better approaches the separable-entangled frontier we are looking for and which is represented in 
Fig. 4. The entangled region in the (a;, y, T)-space is illustrated in Fig. 5 for various values of (a;, a). For T = and 
arbitrary a, the system is entangled if (1 — x)/3 < y < 1 — x, which immediately recovers the Peres' 1/3 criterion for 
x = 0. As T increases, the entangled domain shrinks against the border x + y = 1 and dissapears at a critical T, T c 
(T c = l/ln2 for x = 0). 

If the conditional entropy S q (A\B) is a monotonically decreasing function of q, then the strongest case corresponds 
to the nonnegativity of S 00 (A\B). Whenever the frontier is fully defined by q — > oo, the system is entangled in the 
region defined by 



where 



y* < y < 1 — x and x < x* (15) 



(16) 



y ~ Z-2e- 1 / T {2 + e- 1 / T ~2e- 2 / T ) 
and 

, _ (2(1 - e~ 2 / T )(2 - 3 e~V T + e^' T - e^) - l) y + 1 
X = (5-4 e -2«/T) 



(17) 



However, for sufficiently small T, there is a curved line (corresponding to the nonlinear portion of the frontier in 
Fig. 4), given by y = f(x), which is defined from the conditional entropy for finite q m in and that can be numerically 
determined. In this case, the entangled region is given by 

V* < V < 1 - x and y> f(x). (18) 

An increase of the noise parameter T (see Figs. 4 and 5), i.e., an increase of the relative weight of the high energy 
states of the oscillators, causes shrinking of the entangled region. This is intuitively expected. Also, interestingly 
enough, decreasing the parameter a (see Figs. 4 and 5), which corresponds to an increase of the noise associated with 
the diagonal states, causes the entangled region to enlarge. This is not surprising after all, since this noise makes the 
diagonal terms to become less effective than the off-diagonal ones, which makes quantum coherence to be relatively 
stronger. The interplay of these effects may cause an interesting re-entrance, as illustrated in Fig. 5(a) for the value 
x = 0.6. In this situation, corresponding in fact to small values of a, it is possible to loose the quantum entanglement 
as T increases, and then recover it back, until definitive loss for very large T. To conclude, we believe that the present 
results provide quantitative insight on the influence of the external world (here represented by two different types of 
noise) on a quantum device. 
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FIGURE CAPTIONS 

Fig. 1. S q (A\B) as a function of (x,y) for a = 1, T = 0.1 and typical values of q. 

Fig. 2. S q (A\B) as a function of (x, y) for a — 1, typical values of T and q = 5. 

Fig. 3. The line S q (A\B) — 0, in (y, g)-space, for particular values of x, T — 0.5 and a = 1. 

Fig. 4. Frontier in the (x,y)-space for typical values of T and a— 0.1 (a), 1 (b), 5 (c). The dotted lines correspond 
to the locus of the vertices of the frontier for any T. 

Fig. 5. Frontier in the (x, y, T)-space for fixed values of x and a— 0.1 (a), 1 (b). The entire plane x + y = 1 (VT) 
belongs to the separable-entangled frontier. 
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(a) a = 0.1 
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